Variations in the axisymmetric transport of magnetic elements on 

the Sun: 1996-2010 



David H. Hathaway 
NASA Marshall Space Flight Center, Huntsville, AL 35812 USA 
david . hathawayOnasa . gov 



Lisa Rightmire 

Department of Physics, The University of Alabama in Huntsville, Huntsville, AL 35899 USA 

lar0009@uali.edu 



ABSTRACT 

We measure the axisymmetric transport of magnetic flux on the Sun by cross-correlating 
narrow strips of data from line-of-sight magnetograms obtained at a 96-minute cadence by the 
MDI instrument on the ES A/NASA SOHO spacecraft and then averaging the flow measurements 
over each synodic rotation of the Sun. Our measurements indicate that the axisymmetric flows 
vary systematically over the solar cycle. The differential rotation is weaker at maximum than at 
minimum. The meridional flow is faster at minimum and slower at maximum. The meridional 
flow speed on the approach to the Cycle 23/24 minimum was substantially faster than it was 
at the Cycle 22/23 minimum. The average latitudinal profile is largely a simple sinusoid that 
extends to the poles and peaks at about 35° latitude. As the cycle progresses a pattern of in- 
flows toward the sunspot zones develops and moves equatorward in step with the sunspot zones. 
These in-flows are accompanied by the torsional oscillations. This association is consistent with 
the effects of the Coriolis force acting on the in-flows. The equatorward motions associated with 
these in-flows are identified as the source of the decrease in net poleward flow at cycle maxima. 
We also find polar counter-cells (equatorward flow at high latitudes) in the south from 1996 to 
2000 and in the north from 2002 to 2010. We show that these measurements of the flows are not 
affected by the non-axisymmetric diffusive motions produced by supergranulation. 



Subject headings: Sun: rotation, Sun: surface magnetism, Sun: dynamo 



1. INTRODUCTION 

The structure and evolution of the magnetic 
field in the Sun's photosphere is believed to be 
produced by dynam o processes within the Sun 
(<Charbonneau1l2005l) . This structure and evolu- 



tion must be faithfully reproduced in any viable 
dynamo model. Flux Transport Dynamo (FTD) 
models have recently been used to predict the 
strength of the next solar cycle ( Dikpati et al. I 
2006t lChoudhuri et al. |[2007l ). In these FTD mod- 
els the Sun's axisymmetric flows (differential ro- 



taton and meridional flow) play key roles. The 
meridional circulation transports magnetic flux 
at the surface to the poles, builds up the polar 
fields, and sets the 11-year length of the solar cy- 
cle by its presumed slow equatorward return at 
the base of the convection zone. The differential 
rotation shears the poloidal magnetic field to pro- 
duce strong toroidal fields that erupt through the 
photosphere in sunspots and active regions. 

The structure and evolution of the photospheric 
magnetic field also serves as the inner bound- 
ary condition for all of space weather - condi- 
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tions on the Sun and in the space environment 
that can influence the performance and reliabil- 
ity of space-borne and ground-based technological 
systems. Surface Flux T ransport (SFT) mod els 



have been used since 1984 (jDeVore et al. Ill984h to 



evolve the surface field using the flux that erupts 
in active regions as a source term. This active re- 
gion magnetic flux is then transported across the 
surface by meridional flow, differential rotation, 
and diffusion by supergranules - nonaxisymmet- 
ric, cellular flows that evolve on a time scale of 
about 1-day. The magnetic field structure pro- 
duced in SFT models has been used to model so- 
lar wind structures (wind speed and interplane- 
tary magnetic field) for space weather forecasts 
( Arge fc Pizzo 2009) and to estimate the S un's 
total irradiance since 1713 ( W ang et al. 20051) for 
Sun-Climate studies. 

The strength, structure, and evolution of the 
meridional flow in particular is critically impor- 
tant in both FTD and SFT models. Unfortu- 
nately, the meridional flow is difficult to measure 
due to its weakness. Supergranules have typical 
flow speeds of about 300 m s _1 and differential ro- 
tation has a typical velocity range of ~ 200 m s _1 . 
Yet, the axisymmetric meridional flow has a top 
speed of only 10-20 m s _1 . 

The axisymmetric flows have been measured us- 
ing a variety of techniques. Feature tracking is 
amongst the simplest and oldest but gives differ- 
ent results depending on the nature of the fea- 
tures themselves. Direct Doppler measurements 
can give the plasma flow velocity in the photo- 
sphere but these measurements are subject to sys- 
tematic errors introduced by other solar processes 
and only provide the line-of-sight velocity - which, 
for the meridional flow, vanishes near the equator 
and limb. Global helioseismology provides mea- 
surements of the differential rotation as a function 
of latitude, radius, and time. Local helioseismol- 
ogy can provide measurements of the meridional 
flow as a function of latitude, depth, and time us- 
ing the methods of ring diagram analysis or time- 
distance analysis. 

Sunspots and sunspot groups were amongst 
the earliest featur es used to measu re the ax- 
isymmetric flows. Carrington ( 1859t ) measured 



Newton fe Nunn I (|195lh measured the locations 
of recurrent sunspots groups on successive rota- 
tions as well as individual sunspots on consec- 
utive d ays and found slightl y different rotation 
profiles. iHoward et al. made detailed mea- 

surements of individual sunspot positions recorded 
on photographic plates at Mount Wilson Obser- 
vatory from 1921 to 1982. They found differential 
rotation with uj = 14.52 — 2.84 sin 2 A deg day -1 
(where A is the heliographic latitude) but noted 
that sunspot groups rotate more slowly than in- 
dividual sunspots and large sunspots rotate more 
slowly than small sunspots. 

Sunspots and sunspot group s can also be used 
to measure the meridional flow. iTuominen I ([1942) 
used the latitudinal positions of recurrent sunspot 
groups and found equatorward flow of ~ 1 m s _1 
below ~ 20° latitude and pole ward flow of s imilar 
strength at higher latitudes. IWard (1973) used 
daily sunspot group positions to argue that there 
was no meridional flow at the 1 m s _1 level. How- 
ever, Howard fc Gilman measured the lat- 
itudinal drift of individual sunspots and found an 
equatorward flow of about 3ms" 1 equatorward 
of ~ 25° with an even weaker poleward flow at 
higher latitudes. An obvious drawback to track- 
ing sunspots to measure the axisymmetric flows is 
the limited latitudinal coverage (latitudes <~ 30°) 
and the complete lack of coverage at times near 
sunspot cycle minima. 

Smaller magnetic features, although often con- 
centrated in the active latitudes, do cover the 
entire solar surface and are present even at 



Komm et al~l (|l993Ah 



sunspot cycle minima, 
masked out the active regions in high-resolution 
magnetograms (2048 x 2048 pixel full-disk ar- 
rays) and cross-correlated the remaining mag- 
netic features with those seen the next day from 
1975 to 1991 for several hundred magnetogram 
pairs. They found differential rotation with 
uj = 14.43 -1.77 sin 2 A -2.58 sin 4 A deg day" 1 and 
noted that latitudinal profile was flat ter at sunspot 



cycle m aximum than at minimum. iKomm et al. 



the positions of sunspots on consecutive days 
and noted the presence of an equatorial pro- 
grade current and higher latitude retrograde flow. 



(|1993B( ) used the same technique to measure the 
meridional flow and found a poleward flow that 
varied with sinusoidally latitude, reaching a peak 
velocity of ~ 13 m s _1 at 39° latitude. Fur- 
thermore, they found that the flow speed was 
slower at the sunspo t cycl e maximum than at 
minimum. iMeunuer ( 19991 ) employed this tech- 
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nique (without masking the active regions) us- 
i ng magnetogram pai rs from the MDI instrument 
(jScherrer et al. Ill995t ) on the ESA/NASA SOHO 
mission over the rising phase of sunspot cycle 23 
from 1996 to 1998. She found that the poleward 
meridional flow slowed in the presence of active re- 
gions. In a recent paper lHathawav fc Rightmire 



(2010) did a similar analysis (with masking of the 
active regions) of MDI magnetograms over the 
time period from 1996 to 2009. They obtained 
measurements from over 60,000 image pairs sepa- 
rated by 8-hours. They also found that the merid- 
ional flow was poleward (with a peak velocity of 
r~> llms^ 1 at ~ 45° latitude) and was fast at cycle 
minimum but slow at cycle maximum. In addi- 
tion they noted that the speed of the meridional 
flow was substantially faster at the Cycle 23/24 
minimum than at the Cycle 22/23 minimum. 

Larger magnetic features, and associated struc- 
tures, yield substa ntially different resu l ts for 
the meridional flow. ISnodgrass fe Dailev I (Il996h 
cross-correlated Mt. Wilson coarse array magne- 
tograms (34 x 34 pixel full-disk arrays) obtained 
24-38 days (a solar rotation) apart and found pole- 
ward flow from 10° to 60° but equatorward flow 
at lower latitudes. Their measurements extended 
from 1968 to 1992 - covering three sunspot cycle 
maxima and two minima. They also found a sys- 
tematic dependence of the meridional flow pattern 
on the phase of the solar cycle. Out-flows from 
the sunspot zones were observed to move toward 
the equator in step with the equator ward move- 
ment of the sunspot zones themselves. lLatushkol 
(| 19941 ) used the same low resolution data (after it 
was processed to construct synoptic maps for each 
solar rotation) and also found out-f lows from the 
sunspot zones. ISvanda et al. ( 2007 ) used a mag- 
netic butterfly diagram constructed from synoptic 
maps of the magnetic field averaged over longi- 
tude for 180 equispaced zones in sine-latitude. 
They measured the slope - change in latitude vs. 
change in time - of the magnetic features and 
found a meridional flow with peak velocities of 
about 20 m s _1 at the poleward limit (~ 45°) of 
their measurements. 

Here we measure the axisymmetric motions of 
the the small magnetic elements using the SOHO 
MDI data in which these elements are well re- 
solved. These magnetic elements are presisely 
those whose transport is modeled in SFT models 



and in the surface transport of the FTD models. 

2. DATA PREPARATION 

High resolution full-disk images of the line-of- 
sight magnetic field have been obtained at a 96- 
minute cadence since May 1996 by the SOHO 
MDI instrument. The s e ima ges were used in 
Hathaway fc Rightmire (2010) to find the vari- 
ation in meridional flow strength over solar cy- 
cle 23. They noted in that paper that the MDI 
imaging system appears to be rotated by ~ 0.21° 
counterclockwise with respect to the accepted po- 
sition angle of the Sun's rotation axis. Further- 
more, they found that the accepted position of the 
Sun's rotation axis is in error by ~ 0.08° as was 



noted previously bv lHoward et al. I ([19861 ) and by 



Beck fc Giles I (|2005T ) ■ This small error introduced 
annual variations in the apparent cross-equatorial 
meridional flow. Here we account for those po- 
sitional errors in mapping the full-disk magne- 
tograms to heliographic coordinates by using mod- 
ified values for the position angle and tilt of the 
Sun's rotation axis. In addition, while reprocess- 
ing the data we found a significant reduction in the 
scatter of the measurements if we took the MDI 
image origin to be at the bottom left corner of 
the bottom left pixel - not the center of the pixel 
as indicated in the MDI documentation. Here 



we rep eat the analyses in lHathawav fc Rightmire 



we rep 

4201(1 ) using these corrected magnetic maps and 
examine the variations in both the strength and 
structure of the axisymmetric flows. 

Each full-disk magnetogram is mapped onto he- 
liographic coordinates using bi-cubic interpolation 
onto a grid with 2048 by 1024 equispaced points 
in longitude and latitude for the entire surface of 
the Sun. This mapping gives a close match to 
the spatial resolution of the MDI instrument and 
makes longitudinal and latitudinal velocities lin- 
ear functions of the displacements in the mapped 
coordinates. The line-of-sight magnetic field is as- 
sumed to be largely radial so we divide the mag- 
netic field strength at each image pixel by the co- 
sine of the heliographic angle from disk center to 
minimize the apparent variations in field strength 
with longitude from the central meridian. The 
magnetic fields in sunspots are intense enough to 
produce magnetic pressures similar to the plasma 
pressure (plasma (3 ~ 1). These intense mag- 
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Fig. 1. — MDI magnetogram from 2001 June 
5 04:48 UT mapped to heliographic coordinates. 
Positive magnetic polarities are yellow, negative 
magnetic polarities are blue, and masked areas are 
red. Tickmarks around the border are at 15° inter- 
vals in latitude and in longitude from the central 
meridian. 



netic field elements resist the near-surface plasma 
flow and have their own peculiar motions in lon- 
gitude and latitude which vary depending on the 
size of the sunspot a nd age of the active region 



( Howard et al. 19861 ). For this reason sunspots 



and their immediate surroundings are masked out. 
We found that this could be done quite effec- 
tively by identifying all mapped pixels with field 
strengths \B\ > 500 G and all pixels within 5 
mapped pixels of those points with \B\ > 100 G as 
masked pixels. An example of one of these mapped 
and masked magnetograms is shown in Fig. 1. 

3. ANALYSIS PROCEDURES 

The axisymmetric motions - differential rota- 
tion and meridional flow - of the magnetic ele- 
ments were determined by cross-correlating strips 
of pixels from pairs of mapped images separated 
by 8 hours and finding the shift in longitude and 
latitude that gave the strongest correlation. (Re- 
sults obtained with image pairs separated by 4.8 
hours were substantially the same.) Each strip was 



11 pixels (~ 2°) high in latitude and 600 pixels (~ 
105°) long in longitude. The shift in longitude and 
latitude producing the strongest correlation was 
calculated to a fraction of a pixel by fitting parabo- 
las in longitude and latitude through the correla- 
tion coefficient peaks. This process was performed 
at 860 latitude positions from 75°S to 75°N for 
typically about 400 image pairs over each 27-day 
rotation of the Sun. In all we obtained measure- 
ments from over 60,000 magnetogram pairs. 

The average and the standard deviation of the 
differential rotation and meridional flow velocities 
were calculated at each latitude for each solar ro- 
tation of 27.25 days. The differential rotation and 
meridional flow profiles for each rotation were fit 
with fourth order polynomials in sin A, where A is 
the heliographic latitude. Errors in the fit coeffi- 
cients were estimated using a Monte Carlo method 
with random variations at each latitude character- 
ized by the standard deviations from the measure- 
ments. These polynomial coefficients were also re- 
cast in terms of associated Legendre polynomials 
of the first order. The Legendre polynomial coef- 
ficients are better suited for studies of time varia- 
tions based on the orthogonality of the polynomi- 
als themselves ( Snodgrassl 1984 ). 

The latitudinal profiles of differential rotation 
and meridional flow as measured with these data 
and this method represent the actual axisymmet- 
ric motions of the magnetic elements. Since the 
magnetic elements are fully resolved in these data 
the effects of supergranule diffusion are seen as 
random motions of the magnetic elements and 
these random motions do not introduce any sys- 
tematic errors in our measurements as will be 
shown in Section 7. Profiles were obtained for 178 
rotations of the Sun from June 1996 to September 
2010 with a gap from June 1998 to February 1999 
when radio contact with SOHO was lost and not 
fully recovered. 

4. AVERAGE FLOW PROFILES 

The average differential rotation profile from 
the entire dataset is shown in Fig. 2. The ve- 
locities are taken relative to the Carrington frame 
of reference which has a sidereal rotation rate of 
14.184 deg day -1 . The average differential rota- 
tion profile is well represented by just the three 
terms with symmetry across the equator - 
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V(j>{\) = (a + 6 sin 2 A + csin 4 A) cos A (1) 



with 



35.6 ±0.1 ms~ 



-208.6 ± 1.1 m s" 



(2) 



(3) 



c = -420.6 ± 1.6 m s _1 (4) 
This gives an angular rotation rate profile with 

w(A) = A + Bsin 2 A + Csin 4 A 

with 



(5) 



A = 14.437 ± 0.001 deg day -1 (6) 



B = -1.48 ±0.01 deg day" 1 (7) 



C = -2.99 ± 0.01 deg day -1 (8) 

where coefficient A includes the Carrington rota- 
tion rate. 

This angul ar rotation rate is nearl y identical to 
that found bv lKomm et al. ] (ll993Ah for the time 
interval 1975-1991 using similar data and meth- 
ods. We do find a slight north-south asymmetry 
as seen in Fig. 2 by the deviation of the measured 
profile from the symmetric profile given by the 
dashed line. The differential rotation was slightly 
weaker in the south than in the north. We also 
note the flattening of the profile at the equator 
with a slight (~ 1 m s _1 ) but significant dip from 
±5° to the equator. A similar "dimple" at the 
eq uator was seen previou sly in direct Doppler data 
bv lHoward et al. I (Il980l) and in magnetic element 
motions by ISnodgrass I (|l983l ). 

The average meridional flow profile for the en- 
tire dataset is shown in Fig. 3. Although the aver- 
age meridional flow profile does display substantial 
north-south asymmetry, the profile is well repre- 
sented with just the two anti-symmetric terms - 
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Fig. 2. — The average differential rotation profile 
with the 2a error range for the time interval 1996- 
2010. The symmetric profile given by Eqns. 1-4 is 
shown with the dashed line. 



29.7 ±0.3 m s" 



-17.7±0.7 m s" 



(10) 



(11) 



This gives a peak poleward meridional flow veloc- 
ity of 11.2 m s" 1 at a latitude of 35.2°. This is 
so mewhat slower t han the meridional flow found 
by iKomm et al. (1993B) for the time interval 
1975 to 1991 but with a peak at nearly the same 
latitude. Our average meridional flow profile 
shows substantially different flows in the north 
and in the south. The flow velocity is faster in the 
south and peaks at a higher latitude than in the 
north. The flow in the north appears to nearly 
vanish at the extreme northern limit (75°) of our 
measurements while the flow in the south is still 
poleward with a speed of about 5 m s _1 at the 
southern limit. 

5. VARIATIONS IN FLOW SPEED 

Variations in the amplitudes of the axisymmet- 
ric flow components were examined by plotting 
the rotation-by-rotation histories of the Legendre 
polynomial coefficients. The Legendre polynomi- 
als were normalized so that their maximum values 
were either 1.0 or -1.0. The coefficients that mul- 
tiply them then give the peak velocity for that 
component. The normalized polynomials we used 
are 



V\(X) — (dsin A + e sin 3 A) cos A 



(9) 



pi(A) = cos A 



(12) 



with 
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Fig. 3. — The average meridional flow profile with 
2<T error range for the time interval 1996-2010. 
The anti-symmetric profile given by Eqns. 9-11 
is shown with the dashed line. This profile shows 
substantially different flow in the north and south. 



P 2 X (A) = 2 sin A cos A 



(13) 



/ 1 35 

^3 W — Y 255 (5 sin 2 A — 1) cos A (14) 



P 4 X (A) = 0.947(7sin 3 A-3sinA)cosA (15) 



P 5 X (A) = 0.583(21 sin 4 A -14 sin 2 A + l)cosA (16) 

The Legendre coefficient histories for the differ- 
ential rotation are shown in Fig. 4 along with 
the smoothed sunspot number for reference to 
the phase of the sunspot cycle. The three sym- 
metric components (P/, P 3 \ and P 1 ) dominate 
so we only show the three associated coefficient 
histories. These three coefficients show only a 
slight variation over the sunspot cycle with the 
amplitudes being smaller (less negative - weaker 
differential rotation) at sunspot cycle maximum 
(~ 2002). This "more rigid" differential rotation 
at sunspot cycle maximu m was previously noted 
bv lKomm et al. I (jl993Al ). 

The Legendre coefficient histories for the merid- 
ional flow are shown in Fig. 5 along with 
the smoothed sunspot number. The two anti- 
symmetric components [P^, and P 1 ) dominate so 
we only show the two associated coefficient his- 
tories. These two coefficients show substantial 
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Fig. 4. — The differential rotation associated Leg- 
endre polynomial coefficients (with 2a error bars) 
for the time interval 1996-2010. The coefficient TO 
multiplies P 1 , the polynomial of zeroth order in 
sin A. The coefficient T2 multiplies P3 1 , the poly- 
nomial of second order in sin A. The coefficient T4 
multiplies , the polynomial of fourth order in 
sin A. The smoothed sunspot number divided by 
4 is shown in red for reference. The differential ro- 
tation is slightly weaker (flatter) at sunspot cycle 
maximum. 



variations over the sunspot cycle with the ampli- 
tudes being smaller a t sunspot cycle maximum. 



Ko 



et al. 



(|1993BI ) found similar behavior for 
the time period 1978-1990. 

In addition to this systematic trend over the 
sunspot cycle (fast at minimum and slow at max- 
imum) we find a secular variation in which the 
meridional flow speed was substantially (~ 20%) 
faster at the Cycle 23/24 minimum in 2008 than 
at the Cycle 22/23 mi n imum in 1996. As in 
Hathaway fc Rightmire I ( 2010T ) we note that the 
meridional flow speed was faster for the entire 
interval from 2004 on, than it was at the cycle 
minimum in 1996. This increase in meridional 
flow speed would explain the weak polar fields 
that were produced during that tim e peri od in 



the SFT models o f ISchriiver fc Liu 
(l2009h . 



Wang et al. 



(2008) and 



6. VARIATIONS IN STRUCTURE 

The variations in flow speed shown in the last 
section are produced by and accompanied by vari- 
ations in flow structure. Our analyses produce lat- 
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Fig. 5. — The meridional flow Legendre polyno- 
mial coefficients (with 2a error bars) for the time 
interval 1996-20i0. The coefficient SI multiplies 
the polynomial of first order in sin A. The co- 
efficient S3 multiplies P4, the polynomial of third 
order in sin A. The smoothed sunspot number di- 
vided by 20 is shown in red. The meridional flow 
is slower at sunspot cycle maximum but was even 
faster at Cycle 23/24 minimum in 2008 than at 
Cycle 22/23 minimum in 1996. 



itudinal profiles of the differential rotation and the 
meridional flow for each individual solar rotation 
from June 1996 to September 2010. These profiles 
were obtained at 860 latitude positions between 
±75°. For further analysis we smoothed these pro- 
files with a tapered Gaussian having a FWHM of 
6 latitude points (~ 1°), resampled at intervals of 
1° in latitude, and produced images of these lati- 
tudinally smoothed profiles and of the differences 
between each such profile and the average sym- 
metrized profiles. Little, if any, variation can be 
seen in the full differential rotation profile history. 
However, the meridional flow profile history shows 
substantial variation as shown in Fig. 6. 

The structure of the meridional flow changes 
substantially over the time period represented in 
Fig. 6. The weakening of the poleward meridional 
flow at sunspot cycle maximum (1999-2003) is ev- 
ident in the muted colors surrounding the sunspot 
zones. The strengthing of the meridional flow on 
the approach to Cycle 23/24 minimum in late 2008 
is evident in the intensified colors at most latitudes 
after 2004. 

Fig. 6 also reveals the existence of counter-cells 
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Fig. 6. — The meridional flow profiles for individ- 
ual solar rotations from 1996-2010. Poleward flow 
is indicated by shades of yellow. Equatorward flow 
is indicated by shades of blue. The latitudinal 
centroid of the sunspot area in each hemisphere 
for each rotation is shown in red. The weaken- 
ing of the meridional flow in the active latitudes 
near sunspot cycle maximum is evident as are po- 
lar counter-cells (equatorward flow) in the south 
from 1996 to 2000 and in the north from 2002 to 
2010. 

(equatorward flow). One is found in the south ex- 
tending equatorward to about 60°S at the start 
of the dataset in May of 1996 but that boundary 
moves poleward of our 75° limit by mid-2000. A 
similar counter-cell is seen forming in the north in 
2002 as it dips below 75°N and remains in evidence 
to the end of the dataset in 2010. This long-lasting 
northern counter-cell is clearly the primary source 
of the north-south asymmetry seen in the aver- 
age meridional flow profile (Fig. 3) and may be 
associated with the asymmetry in the differential 
rotation (Fig. 2). The fact that it maintains its 
existence for more than half of the time available 
in this dataset leaves its imprint on the average 
meridional flow profile in the form of the rapid 
drop in poleward flow in the north to near zero at 
75°N latitude. 

Additional details concerning the structural 
changes in the axisymmetric flows are seen when 
the average symmetric flow profiles are subtracted 
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Fig. 7. — The differences between the meridional 
flow profiles for individual solar rotations and the 
average, symmetric profile from 1996-2010. Pole- 
ward flow (relative to the average profile) is indi- 
cated by shades of yellow. Equatorward flow is 
indicated by shades of blue. The latitudinal cen- 
troid of the sunspot area in each hemisphere for 
each rotation is shown in red. The system of in- 
flows toward the sunspot zones is evident as pole- 
ward flow on the equatorward sides of the sunspot 
zones and equatorward flow on the poleward sides. 



from the profiles for each individual rotation. 
These differences from the average for the merid- 
ional flow are shown in Fig. 7. The two counter- 
cells are more obvious here. In addition, these 
difference profiles show a system of in-flows (rel- 
ative to the average meridional flow) toward the 
sunspot zones with poleward (yellow) flows on the 
equatorward sides and equatorward (blue) flows 
on the poleward sides. This suggests that the 
slowdown in the poleward meridional flow seen at 
sunspot cycle maxima is produced by the growing 
strength and latitudinal extent of these in-flows. 

The presence of th ese in-flows was none t heless 
somewhat surprising. Snodgrass fc Dailev ( 19961 ) 
found out-flows from the active l atitudes with 



their low-r esolution magnetic data. IChou fc Dai 
(|200ll ) and lBeck et"al~l (j2002l ) also found out-flows 
from the active latitud es using time-distance helio- 
seismology. However, [Gonzalez Hernandez et al. 



what we see in Fig. 7 using ring-diagram he- 
lioseismology and the structural changes seen in 
the magnetic element motions by lMeunuer (1999) 
also support the presence of these in-flows. 
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(2010) found clear evidence for in-flows much like 



Fig. 8. — The differences between the differential 
rotation profiles for individual solar rotations and 
the average, symmetric profile from 1996-2010. 
Faster (prograde relative to the average profile) 
flow is indicated by shades of yellow. Slower (ret- 
rograde) flow is indicated by shades of blue. The 
latitudinal centroid of the sunspot area in each 
hemisphere for each rotation is shown in red. The 
torsional oscillations are evident as faster flow on 
the equatorward sides of the sunspot zones and 
slower flow on the poleward sides. 

The in-flows toward the sunspot zones are ac- 
companied by the torsional oscillations - varia- 
tions in the differential rotation seen as faster ro- 
tation on the equatorward sides of the sunspot 
zones and slower rotation on the poleward sides 
( Howard fc LaBonte 1980). This is shown in Fig. 
8 by the differences in the differential rotation pro- 
files from the average symmetrized differential ro- 
tation profile. (Note that there are instrumental 
artifacts at the highest latitudes as evident by the 
annual variations in flow speed with faster flow 
near the poles in the hemisphere tilted toward the 
observer. These artifacts may be due to an ellip- 
tical distortion o f the M DI image as reported by 
(|2004| ). However, our efforts to 



Korzennik ct al. 



include this distortion with either the angle they 
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reported or the angle given in the MDI documen- 
tation did not improve the results.) 

These variations in the differential rotation arc 
consistent with the effect of the Coriolis force on 
the in-flows and the counter-cells. Material mov- 
ing equatorward from the higher latitudes will 
spin-down and give slower flows on the poleward 
sides of the sunspot zones while material moving 
poleward from the equator will spin-up and give 
faster flows on the eq uatorw a rd sid es. This scene- 
no was sugg ested by ISpruid ([20031) as a response 
to cooling in the sunspot zones by excess ther- 
mal emission from faculae. Earlier, ISnod grass 



(|1987l ) had suggested that in-flows and the tor- 
sional oscillations were part of a system of az- 
imuthal convection-rolls which migrate equator- 
ward during each sunspot cycle. These convection- 
rolls should have out-flows at some undetermined 
depth below the surface - a possible source of the 
out-flows seen in some of the helioseismology stud- 
ies. The Coroilis force acting on the long-lasting 
northern counter-cell should slowdown the rota- 
tion at the affected latitudes. This may be the 
source of the north-south asymmetry in the aver- 
age differential rotation profile (Fig. 2). 

7. EFFECTS OF DIFFUSION ON FLOW 
MEASUREMENTS 

The magnetic elements under study here are 
also subject to a diffusion-like random walk by 
the nonaxisymmetr ic cellular flows - supergran- 
ules in particular ([Leighton 
dom walk transports the weak 



1964|). 



This ran- 
magnetic ele- 
ments in both longitude and latitude and leads 
to the formation of large unipolar areas from 
the precedin g and following magnetic flux in ac- 



tive regions (pmithsonl 119731 ). This random walk 



might contribute to the meridional flow we mea- 
sure due to resultant changes in the magnetic 



pattern. In SFT models ( 


DeVore et al. 


1984: 


van Ballegooiien et al. 1998 


; IWang et al. 


2002. 


2005. 20091 Schriiver & Liu 12008) this process 



is represented by a diffusivity coupled with the 
Laplacian of the magnetic field. We would expect 
that this might produce a meridional flow signal in 
the form of out-flows from the sunspot zones where 
the magnetic field is concentrated. Although what 
we observe is actually in-flows toward the sunspot 
zones, the effects of diffusion might nonetheless al- 



ter the structure and evolution of the meridional 
flow we measure. Given this caveat, we undertook 
an investigation of the effects of supergranule dif- 
fusion on our measurements. 

Hathaway et al. (2010) have recently produced 



a model of the photospheric flows which includes 
the cellular flows, supergranules in particular, ob- 
served with the SOHO MDI instrument. The cel- 
lular flows in this model have velocity spectra, life- 
times, and motions that match those seen in the 
MDI data itself. We have taken the vector veloc- 
ities from this model and used them to transport 
magnetic elements whose initial spatial distribu- 
tion was taken from an MDI synoptic magnetic 
map. We then used our analysis procedures to 
measure the axisymmetric flows. We isolated the 
effects of diffusion by only including the evolving 
cellular flows. We do not include the axisymmet- 
ric meridional flow or differential rotation and the 
cellular flow pattern itself does not participate in 
any axisymmetric meridional flow or differential 
rotation. 

The cellular flow simulation produced vector 
velocities on a heliographic grid with 4096 by 
1500 equispaced points in longitude and latitude 
from an evolving velocity spectrum that extended 
to spherical wavenumbers of 1500 (supergranules 
have spherical wavenumbers of ~ 100). The ini- 
tial magnetic field distribution was taken from an 
MDI synoptic magnetic chart for Carrington rota- 
tion 2000 (mid-2003 - just after the peak of the 
sunspot cycle). Our magnetic flux transport sim- 
ulation was calculated on a grid the same size as 
our mapped magnetograms. At each pixel in our 
simulated magnetic map we introduced a number 
of 1000 G magnetic elements with filling factors of 
5% until the average field strength in that pixel 
equaled the observed field strength (a single cle- 
ment in a pixel would produce a field strength of 
50 G). This process required some 120,000 mag- 
netic elements. These elements were then trans- 
ported explicitly by the velocity field from the cel- 
lular flow simulation in 15-minute time steps for 
10 days. 

Examples from the simulated magnetic maps 
are shown in Fig. 9. The magnetic elements 
are transported to the borders of the cells and 
then continue to move as the cells themselves 
ev olve. (This was show n in previous simulations 
bv ISimon et al. ( 2001 ).) The magnetic elements 
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Fig. 9. — Simulated magnetic map regions at 1- 
day intervals. These regions were extracted from 
the full simulated magnetic maps at the start of 
days 1-5 from an area bordered by the equator, 
60°N, and longitudes 109° and 126°. The evolv- 
ing magnetic network is evident in the changing 
magnetic structures. 

retain their identities throughout the simulation 
and do not interact with each other. If opposite 
polarities occupy a pixel they do cancel each other 
in terms of the mapped magnetic field strength but 
they continue to retain their identities and move 
with the simulated flow. 

These magnetic maps were processed with the 
same analysis procedures used with the MDI mag- 
netic maps by selecting a "central meridian" lon- 
gitude and correlating strips of pixels with those 
from a map 8-hours later. This was done for a 
series of cental meridians at 1-hour intervals over 
the 10 simulated days. This resulted in 559 mea- 
surements of the axisymmetric flows covering the 
full range of longitudes and the full 10 days. Fig. 
10 shows the meridional flow measured from these 
magnetic maps. The results have similar noise lev- 
els to single rotation averages from MDI but show 
no evidence of any systematic meridional flow. 

8. CONCLUSIONS 

We have measured the axisymmetic motions of 
magnetic elements on the Sun by cross-correlating 
strips of data from magnetic maps acquired at 
96-minute cadence by the MDI instrument on 
SOHO. Our measurements cover each rotation 
of the Sun from June 1996 to September 2010 
with the exception 8 rotations when the data 






Latitude 

Fig. 10. — Meridional flow profile measured from 
magnetic features subjected to random walk by 
non-axisymmetric cellular flows. Our meridional 
flow measurements do not include any systematic 
errors due to these random (and spatially resolved) 
motions. 



were unavailable. Although we exclude the 
magnetic elements in sunspots themselves, the 
magnetic elements we track are in fact those 



lar fields in SFT models ( 


DeVore et al. 


1984; 


van Balleeooiien et al. 1998 


; Wang et al. 


2002. 


2005. 200*} Schrijver& 


Au 


2008) and in FTD 


models (iDikpati et al. 1 


2006: Choudhuri et al. 



120071 ). With these data these magnetic elements 
are well resolved and the random motions due to 
supergranules appear as just that - random mo- 
tions that do not alter our measurements of the 
axisymmetric flows. 

The differential rotation we measure agrees well 
with previous measurements using similar data 



and methods (jKomm et al. I [1993AJ) . Although 



the average differential rotation profile is slightly 
asymmetric this asymmetry may be specific to 
the time period and the presence of the merid- 
ional flow counter-cell in the north. The tor- 
sional oscillation signal (Fig. 8) compares well 
with the near surfa ce pattern from helioseismology 
( Howe et al. 1 12009) and does not require averaging 
the two hemispheres together. 

The meridional flow we measure also agrees well 
with previous measurements using s i milar data 



and methods (Komm et al. 1993B: Mcunuer 



1999h but with interesting differences and more 
detail. The average meridional flow speed we 
found from 19 96 to 2010 was somew hat slower 
than found by iKomm et al~l (|l993Bl ) from 1978 
to 1991. We both find that the flow is faster at 
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cycle minima and slower at maxima. Here we find 
that this slow-down can be attributed to a system 
of in-flows toward the sunspot zones which, when 
superimposed on the average meridional flow pro- 
file, lowers the p eak flow velocity at cycle maxima 



( Meunuer I ll999ft . Our slower average meridional 



flow speed is somewhat surprising since our data 
included two (fast) minima and one maximum 



while the iKomm et al. (1993B) data included 



two (slow) maxima and one minimum. 

An important difference for understanding the 
long, drawn-out, and low Cycle 23/24 minimum 
is the faster meridional flow after 2004 compared 
to the flow at the Cycle 22/23 minimum in 1996. 
This faster meridional flow produces weake r polar 
field s in the SFT models of lSchriiver fc Liu I (|2008h 



and IWang et al. (2009). Weaker polar fields pro- 



duce weak followi ng cycles which typically have 



long, low minima (jHathawav Il2010h 



In spite of this agreement, our average merid- 
ional flow profile is problematic for the SFT mod- 
els. All of the SFT modeling groups use merid- 
ional flow profiles which peak at low latitudes 
or do not extend poleward of 75°. Comparisons 
between our symmetrized profile and those used 



in th r ee SFT calculation s (Ivan B allcgooiicn et al. 



19981: IWang et al. 1 120091: ISchriiver fc Title 1 12001) 



are shown in Fig 11. 
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Fig. 11. — Symmetrized meridional flow profile 
from this paper (solid line) plotted with merid- 
ional flow profi les use in the Surfa ce Flux Trans- 
port models of IWang et al" ]j2009h (dashed line) 



van Ballcgooii en et al. I (|1998|) (dotted line) and 



Schriiver fc Title I (|200lh (dot-dashed line). 



All three SFT profiles fall below our measured 
profile at th e high er latitudes - above 30° for 
Wang et al. I (l2009l) . 45° for ISchriiver fc Title 



(|200l[ ). and 60° for lvan Ballegooiien et aU (|l998T ). 
Using our average meridional flow profile in these 
models without compensating processes leads to 
polar fields substantially stronger than those ob- 
served. Compensating processes might include the 
counter-cells along with the north-south asymme- 
try or neglected physical proce sses - for exam- 
ple ra dial diffusion suggested bv iBaumann et ah 
(200§). 



The nearly 20% change in meridional flow speed 
from Cycle 22/23 minimum in 1996 to Cycle 23/24 
minimu m in 2008 is problematic for t he FTD 
models. Dikpati fc Charbonneau ( 19991 ) showed 
that with their FTD model increasing the surface 
meridional flow speed from 2 m s^ 1 to 20 m s" 1 
changed the surface polar field strength from 130G 
to 350G while changing the cycle period from 77 
years to 11 years. The faster meridional flow in 
this model should have produced a shorter cycle 
with stronger polar fields. Yet, observations reveal 
a very long cycle with much weaker polar fields. 

We have shown that our data, with its high 
spatial resolution and rapid cadence, fully resolve 
the magnetic element motions produced by super- 
granule "diffusion" and thus yield measurements 
of the meridional flow with out any systematic er - 



rors due to that diffusion. IKomm et al. I (|19 93B) 
used data with similar spatial resolution but lower 
cadence (daily rather than hourly) and found 
similar results. Howe v er, Snodgrass fc Dailev 



(<1996l) and lLatushko I (Il994l) used data with 



much lower spatial resolution and much longer 
time-lags (monthly) and found significant dif- 
ferences. These low spatial resolution data do 
not resolve the individual magnetic elements. 
They image the emsemble magnetic patches 
whose motions do include the effects of diffu- 
sion. We suspect that the magnetic pattern 
diffusion gave the e quatorial flows at low lati- 
tudes measured by ISnodgrass fc Dailev" ( 19961 ) 
and the out-flows from t he su nspot zones seen 
bv ISnodgrass fc Dailev I (|19961 ) and lLatushko 



(119941) . and mor e rapi d high-latitude flow seen 



bv lSvanda et al. I (|2007l ). 

Comparisons of our measurements with those 
from other data types (direct Doppler velocities, 
sunspot motions, and helioseismology) are sub- 
ject to problems associated with the characteristic 
depth of the measurements. The Sun has a sur- 
face shear layer produced largely by the granule 
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and supergranule flo ws which tend to conse rve 



angular momentum ( Foukal k Jokipii Il975l ) 



slowing down the rotation of the surface layers 
and speeding up the rotation down to depths of 
about 35 Mm. This inward increase in rotation 
rate should be accompanied by an inward de- 




creas e in the meridional flow s peed ([H athaway 
■ a feature noted by lHathawav et al. 



198 
(201 



ules. 



in the meridional motion of supergran- 
This is consistent with the slower rota- 



tion rate and faster meridional flow seen in di- 
rect Doppler measurements representative of the 



photosphere (jUlrich et al. Ill988l : lUlrich lEzOloj ) 



as- 



suming that the magnetic elements are rooted 
in somewhat deeper layers. Sunspots should be 
rooted even deeper yet and sunspot s show rota- 
tion rates which are even more rapid ( Ward 19661: 
Howard et al~lll986l ) and meridional motions that 
are v anishingly small (IWard _ 19731) 



are v amsmngiy small II ward I naval or equator - 
ward (jTuominen Ill942t iHoward fc Gilman l ll986) 



While helioseismology studies indicate both out 



flows ( Chou k Dai 2001; Beck et al 



in-flows (Gonzalez Hernandez et al. 



2002) and 
2010h . this 



may be due to differences in both the methods 
used and the associated depths of the measure- 
ments. Helioseismology does provide support- 
ing evidence for the variatio ns in meridional flow 
speed over the sunspot cycle ( Basu k Antia f 2003: 



Gonzalez Hernandez et al. Il2010l ). 

Our observations of in-flows toward the sunspot 
zones may help us understand the origins of the 
torsional oscillations. The strength and structure 
of these in-flows are go od mat c hes t o the flows 
predicted in the model of lSpruitl (|2003l ). However, 
helioseismology indicates that the torsional oscilla- 
tions may originate well below the surface at high 
latitudes ([Basu k Antia 1120031) and thus may not 
be forced by the effects of localized surface cooling. 

Finally, we reitterate our point that the mag- 
netic elements whose motions we study are pre- 
cisely those elements whose transport is modeled 
in SFT models and at the surface in FTD mod- 
els. Both SFT and FTD models must employ the 
measured axisymmetric transport of those mag- 
netic elements to conform with observations. 
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